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Abstract 

This paper is the second part of a study of the quantum free particle on spherical 
and hyperbolic spaces by making use of a curvature-dependent formalism. Here we 
study the analogues, on the three-dimensional spherical and hyperbolic spaces, S% 
(k > 0) and (k < 0), to the standard spherical waves in E 3 . The curvature 
k is considered as a parameter and for any k we show how the radial Schrodinger 
equation can be transformed into a K-dependent Gauss hypergeometric equation 
that can be considered as a K-deformation of the (spherical) Bessel equation. The 
specific properties of the spherical waves in the spherical case are studied with great 
detail. These have a discrete spectrum and their wave functions, which are related 
with families of orthogonal polynomials (both ^-dependent and k- independent), and 
are explicitly obtained. 
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1 Introduction 



This article can be considered as a sequel or continuation of a previous paper [1], which was 
devoted to the study of the quantum free particle on two-dimensional spherical and hyperbolic 
spaces making use of a formalism that considers the curvature n as a parameter. Now, we 
present a similar analysis but introducing two changes related with the dimension of the space 
and with the states of the quantum free particle we are looking for. Now we work in a three- 
dimensional space, and we look for the states analogous to the Euclidean spherical waves, which 
are determined among all free particle states by the condition of being separable in the geodesic 
polar coordinate system. We follow the approach of [1], which contains the fundamental ideas and 
motivations, and we also use the notation, ideas, and results discussed in some related previous 
studies [2]- [5]. 

There are two articles that are considered of great importance in the study of mechanical 
systems in a spherical geometry (see [1] for a more detailed information; we just make here a 
quick survey in a rather telegraphic way). Schrbdinger studied in 1940 the hydrogen atom in 
a spherical space [6] and then other authors studied this problem (hydrogen atom or Kepler 
problem) [7]-[12] or other related questions (as, e.g., the quantum oscillator on curved spaces) 
[13]- [16]. Higgs studied in 1979 the existence of dynamical symmetries in a spherical geometry 
[17] and since then a certain number of authors have considered [18] -[39] the problem of the 
symmetries or some other properties characterizing the Hamiltonian systems on curved spaces 
(the studies of Schrodinger and Higgs were concerned with a spherical geometry but other authors 
applied their ideas to the hyperbolic space). In fact these two problems, the so-called Bertrand 
potentials, have been the two problems mainly studied in curved spaces (at the two levels, 
classical and quantum). Nevertheless in quantum mechanics there are some previous problems 
that are of fundamental importance as, for example, the quantum free particle or the particle in 
a spherical well. 

This article is concerned with the study of the spherical waves for a quantum free particle on 
spherical and hyperbolic spaces (an analogous problem was studied in [1] but in three dimensions 
and making use of K-dependent parallel geodesic coordinates, which affords the analogous of plane 
waves). There are some questions as, for example, (i) analysis of some ^-dependent geometric 
formalisms appropriate for the description of the dynamics on the spaces with constant curvature 
k, (ii) transition from the classical ^-dependent system to the quantum one, (iii) analysis of the 
Schrodinger separability and quantum superintegrability on spaces with curvature that have 
been discussed in [1], and therefore they are now omitted (or revisited in a very sketchy way). 
Thus, this paper is mainly concerned with the exact resolution of the K-dependent Schrodinger 
equation, existence of bound states, and with the associated families of orthogonal polynomials. 

In more detail, the plan of the article is as follows: In Sec. 2 we first study the Lagrangian 
formalism, the existence of Noether symmetries and Noether momenta and then the K-dependent 
Hamiltonian and the quantization via Noether momenta. In Sec. 3, we solve the K-dependent 
Schrbdinger equation and then we analyze with great detail the spherical n > case, writing 
explicitly the spherical waves on a 3D-sphere and discussing their Euclidean limit when the 
curvature of the sphere goes to 0. The study of the hyperbolic k < case is only sketched, but 
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the details displays several important differences which would require a separate study. Finally, 
in Sec. 4 we make some final comments. 



2 Geodesic motion, ^-dependent formalism and quan- 
tization 



We first present a brief introductory comment on some possible approaches to the two-dimensional 
manifolds with constant curvature re: the sphere S 2 (re > 0), Euclidean plane E 2 , and hyperbolic 
plane H 2 (re < 0), and then we move to the corresponding three-dimensional spaces: the sphere 
S 3 . (re > 0), Euclidean space E 3 , and hyperbolic space (re < 0). 

If we make use of the following re-dependent trigonometric (hyperbolic) functions 



cos yre x if re > 0, 

Ck(x) = I 1 if K = 0, 

s cosh^/ — re x if re < 0, 



S K (x) = 



7^ sin y^x if k > 0, 
x if re = 0, 

J= sinhv/— rex if re < 0, 



then the expression of the differential element of distance in geodesic polar coordinates (r, cf)) on 
M 2 = (S* 2 ,E 2 ,F 2 ) can be written as follows 

dl 2 R = dr 2 + S 2 K (r)d<t ) 2 , 

so it reduces to 

dl\ = dr 2 + (sin 2 r) d<f> 2 , dl% = dr 2 + r 2 d<f> 2 , dl 2 _ 1 = dr 2 + (sinh 2 r) d<f> 2 , 

in the three particular cases re = 1, 0, —1 of the unit sphere, the Euclidean plane, and the 'unit' 
Lobachewski plane. If we make use of this formalism then the Lagrangian of the geodesic motion 
(free particle) on M 2 is given by [5, 28, 29] 
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(1) 



JL(re) = (-)( V 2 + S^) . 
Now if we consider the re-dependent change r — > s = S K (r) then the Lagrangian IL(re) becomes 

and, if we change to 'cartesian coordinates for s' defined as x = s cos (j),y = scos^, we arrive to 



L ( K ) = 2 ( l - K s 2 ) v l+ v y- K ( xv y-y v *) 2 



2 2,2 

s z = x z + y. 



that is the Lagrangian studied in Ref. [1, 4, 5] (the relation of this Lagrangian with the Lagrangian 
of Higgs is also discussed in [1, 5]). 
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We notice that in the sphere case, in addition to the usual geodesic polar coordinate singularity 
at the origin (the 'North pole') r = 0, which passes to the coordinate s at s = 0, the chart (s, 4>) 
covers only the upper hemisphere < r < tt/(2^/k) as the coordinate s ceases to be related to r 
on a one-to-one basis at the equator r = 7r/(2y / K), where S K (j) reaches its maximum; however 
the lower hemisphere 7r/(2y / «) < r < ir / 'yfk can be also covered by another chart, with s still 
given by s = S K (f) with a singularity at s = 0, which on the lower hemisphere corresponds to 
r = 7t/a/k, the point antipodal to the origin (the 'South pole'). 



2.1 Lagrangian formalism, Noether symmetries and Noether mo- 
menta 

Let us start with the following expression for the differential element of distance dl K in the 
family M% = (S^,1E 3 ,H 3 ) of three-dimensional spaces with constant curvature k written in 
(s,6,(j)) coordinates (recall s is not the geodesic radial coordinate): 



dlt 



ds 2 



+ s 2 d9 2 + s 2 sin 2 9 dc/) 2 



(2) 



1 — K S 2 

reducing in the particular cases of unit sphere, Euclidean plane, and 'unit' Lobachewski plane to 

ds 2 



dlf 



+ s 2 (d9 2 + sin 2 ## 2 ), 



1+ s 2 

Then the following six vector fields 

Xi(k) = y/l — ks 2 (sin 9 cos 4 
X 2 (k) 



1 - s 2 

dll = ds 2 + s 2 {d9 2 + sin 2 9 # 2 ) , 
dl 2 = 



+ s z (d9 2 + sin 2 9 <t> z ) . 



d 1 r , . d .sine!), d , 

— + - cos 9 cos — - ^ — 
l os s 09 sm9 ocp 

r, o \, ■ n ■ ,\ ° 1 r/ n ■ ,\ & ,cos0. d . 

v 1 k s (sin^ sin 0) — + - (cos0sm<£) — + (^-^-) ^- 
L as s 09 sin 6* a0 

9 1 . „ 



(cos 0) - sin — 

os s o9 



and 



v ■ a d r cos4> d d sm<t> d d 

Yi = - sm *d9~ W^' F2 = cos ^- W^' F3 = ^' 



are Killing vector fields, that is, the infinitesimal generators of isometries of the K-dependent 
metric dl 2 . 

The Lie brackets of the vector fields Xi(n) are given by 

[X 2 (k),*i(k)] = X K Y 3 , [X 3 (k),X 2 (k)} = \ K Y-i , [Xi(k),X 3 (k)] = X K Y 2 , 
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with A K given by A K = k\/\ — k s 2 . The other Lie brackets are K-independent and similar to the 
Lie brackets of the Euclidean case; that is 

[Y 2 ,Y 1 ] = Y 3 , [Y 3 ,Y 2 ] = Y 1 , [Y 1 ,Y 3 ] = Y 2 , 

and so on. All these Killing vector fields close a Lie algebra that is isomorphic to the Lie algebra of 
the group of isometries (either SO(4), I SO (3), SO(l, 3)) of the spherical, Euclidean or hyperbolic 
spaces depending of the sign of k. Notice that only when k = (Euclidean space), the vector 
fields Xi(n), i = 1,2, 3, will commute between themselves. 

Now, let us consider the geodesic motion on M 3 , that is, the dynamics determined by a 
Lagrangian L, which reduces to the K-dependent kinetic term T(k) without a potential 

L = T(k) = (\) + s 2 v 2 + s 2 sin 2 Ov 2 ) , (3) 

where the parameter k can take both positive and negative values. We already mentioned that 
in the spherical case the coordinate chart we are dealing with covers only the 'upper' half- 
sphere; we see that this Lagrangian becomes singular at the 'equator' where r = and hence 

1 — ks 2 = 0, so in this case the study of the dynamics will be restricted to the interior of the 
interval < s < 1/\/~k which corresponds to the upper half sphere. 

The Lagrangian L = T(k) is invariant under the action of the the K-dependent vector fields 
Xi{n) and Yj, i = 1,2,3, in the sense that, if we denote by X\(n) and the natural lift to 
the tangent bundle (phase space TM% with M% representing S%, E 3 , or H%) of the vector fields 
Xi(n) and Y«, i = 1, 2, 3, then the Lie derivatives of T(k) vanish, that is 

Xft«)(T(«)) = 0, 5?(T(k)) = 0, i = 1,2,3. 

They represent six exact Noether symmetries for the geodesic motion. If we denote by 0l the 
Lagrangian one-form 

2 ) ds + s vgdO + s sin Ov^dcf) , 

then the associated Noether constants of the motion are given by the following: 
(P) The three functions P\(k), P 2 (k), and Ps(k), defined as 

p i (K) = i(xf( K ))e L , i = 1,2,3, 

that are ^-dependent and given by 

D I 

Pi(k) = (sin 9 cos (f>) , s = = + (s y 1 — ks 2 ) [(cos cos <p) vq — (sin 9 sin ef>) v^] , 

V 1 — K S 

u 1 

P2{k) = (sin 9 sin (f>) —=^= + (syl — ks 2 ) [(cos 9 sin 4>)vq + (sing cos^>)^] , 

P3(ac) = (cosg)— s — (sy I — k s 2 )sva9 vg . 

V 1 — K S 2 
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(J) The three functions J\, J2, and J3, defined as 

Ji(K) = i(Y?)9 L , i = 1,2,3, 



that are ^-independent functions and given by 

Ji = — s 2 (sin (f> v + sin 9 cos 9 cos <j> v^) , 
J2 = s 2 (cos0t>0 — sin^cos^sin^t;^) , 
J3 = s 2 sin 2 -u^ . 



2.2 ft- dependent Hamiltonian and Quantization 

The standard method for the quantization of a Hamiltonian on a Riemannian manifold is to 
make use of the Laplace-Beltrami operator for the free part (kinetic energy) of the Hamiltonian. 
Nevertheless we recall that the standard procedure in a Euclidean space is to first quantize the 
momenta as self-adjoint operators and then, making use of the quantum momenta, to obtain the 
quantum version of the Hamiltonian. Our idea is to translate this momentum-approach to the 
case of spaces with curvature k but changing the quantization of the canonical momenta by the 
quantization of the Noether momenta which are taken as the basic objects (this is one of the 
reasons why we have studied the properties of the Killing vectors and Noether momenta with 
great detail). So, we present the quantization of the system in two steps: (i) quantization of the 
Noether momenta as self-adjoint operators and then (ii) quantization of the Hamiltonain making 
use of the quantum Noether momenta. 

The Legendre transformation (s, 9, (/>, v s , v$, v^) — > (s, 9, (f),p s ,P9,Pcf>) is given by 



Ps 



1 



Pe 



2 2 ■ 2 a 

s vg, p<t, = s sin 9 v^, , 



so that the expression of the K-dependent Hamiltonian turns out to be 

^2 



v '••-^ • • si ;,2, 



(4) 



The six Noether momenta become 



Pi(k) = \/l — ks 2 (sin 9 cos <p) p s H — [(cos 9 cos < 



,sm f 



s 
1 



)Pe 



sin 1 

, COS 1 



-)pA 



Pi{k) = \A - « s 2 (sin 9 sin <f>) p s + - [(cos 9 sin 4>)pe + {—^~e )P<j>] 

s Sill f / 

P^in) = y/l — ks 2 (cos9)p s sin9pg 



and 



.cos < 



J 1 = - sm<f)p e - (- ^)pd>, J2 = cos (pp g 

tan o 



, sm <p A 
v tan0 y 



■ , J3 — P(j>, 



with Poisson brackets 



{Pi(K),H( K )} = 0, {J 1 ,H(k)} = 0, 



1,2,3, 
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and 

{P 1 (k),P 2 {k)} = kJ 3 , {P 2 (k) , P 3 (k)} = k Ji , {P 3 (k),P 1 (k)} = kJ 2 . 
The other Poisson brackets are similar to the Poisson brackets of the Euclidean case; that is, 

{Pi(«),Ji} = 0, {P 1 (k),J 2 } = P 3 (k), {P 1 (k),J 3 } = -P 2 (k), 

and so on. Note the change of the order in the Poisson brackets. This is motivated because of 
the property [Xf,X g ] = —X^ g y, that is, the map [,]—>■{,} is a Lie algebra isomorphism but 
with a change of the sign. 

Making use of this formalism, the Hamiltonian of the K-dependent oscillator can be rewritten 
as follows 

(5) 
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H(K) = ( — ) P( + P 2 2 + Pi + K (jf + 4 + J^) 



Proposition 1 The only measure that is invariant under the action of the three vector fields 
Xi(n) and the three vector fields Yi, is given in coordinates (s,6,<f)) and up to a constant factor 
by 



( s smO \ 



dOi 



This property is proved as follows. The most general expression for a volume three-form is given 
by 

uj = n(s, 6, 4>) dsAd9Ad(fr 
where /j,(s, 6, <f>) is a differentiable function to be determined. Then the conditions 

Cy i dco = , Cy i dw = , i = 1, 2, 3, 

lead to the following value for the function fi: 

s 2 sin 9 



K( S ™^_ ), 



where K is an arbitrary constant. Assuming K = 1 we obtain dfi K . 

This property suggests the appropriate procedure for obtaining the quantization of the Hamilto- 
nian H(k). The idea is to work with functions and linear operators defined on the space obtained 
by considering the three-dimensional space endowed with the measure dfx K . This means, in the 
first place, that the operators Pi, P 2 , and P 3 , representing the quantum version of of the Noether 
momenta Pi, P 2 , an P3 must be self-adjoint not in the standard space L 3 (IR 3 ) but in the space 
L 2 (IR 3 ,^ K ). If we assume the following correspondence: 

/ ■ n ^9 1 n d .sin <p. d 

(sine cos 0) tt- + - cos 6 cos 0) — - (— - — 
os s 06 smO aq> 

Q 1 d cos <p d 

(sin6»sin0) — + - [(cos sin 0) — + {-^) ttt 
os s 06 s\nB d(f> 



p 1 ^ p 1 = -ih^/T^ 
p 2 ^ % = -ih^/T^ 



K S^ 



K S^ 
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and 



then we have 



Pi + P| + P$i->- 



J l !->■ J I 
J 2 •->■ ^2 
^3 •-»• ^3 

(1 - KS 2 



— i h 

— ih 



d cos N d 
sin — + ( 



5(9 v tan0' 90 
sin d 

""+88- 



a 2 



+ 



90' 
1 / 9 2 



+ 



9 2 1 
+ 
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9s 2 s 2 Vd6> 2 sin 2 0d0 2 tan#<90 



+ 



2-3ks 2 d 
r ds 



and 



J{ + J| + J| i-> - ft 2 



9^ | 1 
dfl 2 + sin 



2 fld0 2 



+ 



1 d 

tan6d6 



in such a way that the quantum Hamiltonian H(n) 



H{K) = (^") + ^ + V + K (Jl 2 + J 2 2 + J3 2 ; 

2m L 



is represented by the following differential operator: 
H(k) 



— In- 2 \— 2 ~ 3ks2 9 1 / d 2 1 



8 2 1 
+ 



(6) 



(7) 



We note that this operator is self-adjoint with respect the measure d\i K and also that it satisfies 
the appropriate Euclidean limit (in this limit s goes to the Euclidean radial coordinate r, so to 
conform with the standard Euclidean usage we write r in this expression): 



lim H(k) 



2m 



d 2 2d 1 / d 2 
+ - 7^ + — ' 



+ 



f 



d 2 1 

+ 



d 



Q r 2 r Q r r 2 \ QQ2 g j n 2 q Q^p, £an Q QQ 



Note that if we write H{n) = #1 + # 2 + #3 with H = Pi +kJ { , i = 1, 2, 3, then [H ,Hj}^0, 
i 7^ j, because of the At-dependent terms. Nevertheless [Hi , Hjk] = 0, Hjk = Hj + H^, i ^ j ^ k, 
so that H(k) can be written as sum of two operators that conmute in several different ways. 
Finally we also note that H(k) can also be written as 

H(k) = H p + kHj, [Hp,Hj]=0, H P = P 2 , H J = J 2 , 

that corresponds to the approach considered in this paper. 
We close this section with the following observations: 
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1. Only for reference, we mention that had we used the polar geodesic coordinates (r,9,q 
then the Hamiltonian would have be represented by the following differential operator: 



H(k) 



2m 



Sl(r)dr 



d 



+ 



i / d 2 i 



+ 



d 2 



+ 



dr J S 2 {r) \d9 2 sin 2 6d(t) 2 tan 9 89 



(8) 



2. The measure d/j, K was introduced as the unique measure (up to a multiplicative constant) 
invariant under the Killing vectors. We have verified that it coincides with the correspond- 
ing Riemann volume in a space with curvature k. 

3. The Noether momenta quantization procedure was motivated in the first paragraph of 
section 2.2. At this point we must clearly state that the final result (that is, the expression 
of the Hamiltonian operator) coincides with the one that would be obtained making use of 
the Laplace-Beltrami quantization. In fact, the fundamental point for the validity of our 
approach was that the Hamiltonian is the quadratic Casimir of the isometry algebra; this 
leads to the Laplace-Beltrami operator. 



3 ^-dependent Schrodinger equation 

The Schrodinger equation: 

H(k)^ = E^, E = E p + kEj, (9) 
leads in the coordinates (s, 9, 4>) we are using to the following ^-dependent differential equation: 



2m 



(I- us 1 



2 



+ 



2-3ks 2 d 



+ 



1 / d 2 Id 2 Id 



+ 



+ 



* = E$> 



(10) 



ds 2 ' s ds s 2 \d6 2 sin 2 6 dcj) 2 tantfdfl, 
Thus, if we assume that ^(s, 9, 4>) can be factorized in the form 

*( 8 ,0,<l>) = R{8)Y Lm (0,<l>), (11) 
where R is a function of s and Yi, m {9, 4>) are the standard K-independent spherical harmonics 



( d 2 1 d 2 1 d\ 



\89 2 sui 2 9d<p 2 t<m9d9 
then we arrive to the following K-dependent radial equation: 

2-3ks 2 d L(L + 1 



h2 \m 2^ d 2 
'2m-[ {l - Ks) di 2 + 



ds 



L(L + l)Y Lr , 



R = ER, R = R(s) 



that can be rewritten in the form 



(1 - ks 2 



ds 2 



2-3ks 2 d L(L + 1) g2 
s ds s 2 



R = 0, 



(12) 
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where £ is defined for any value of the curvature so that it bears with the energy the same 
relation as the modulus of the wave vector k has with E in the Euclidean case: 

2 _ 2mE 

If we adimensionalyze the radial variable s and the curvature k to the new adimensional variables 

c ~ /c2 2 ~ 2 

p = ts, k = k/c , ks =np , 

then we arrive to the following equation for the function R(p): 

P 2 {1 -np 2 ) R" + p(2 - 3k p 2 ) R! + (p 2 - L(L + 1)) R = 0, (13) 

that represents a K-dependent deformation of the spherical Bessel differential equation 

p 2 R" + 2pR' + [p 2 — L(L + 1)]R = 0. (14) 

This 'deformed' spherical Bessel equation can be solved in power series using the method of 
Frobenius. First the function R must be written as follows 

R = p»f(p,K), 

and then it is proved that p must take one of the two values p± = L or p2 = —L — 1. Choosing 
p = L, in order to translate to the unknown function f(p) the condition R has to satisfy to be 
well defined at the origin in a simpler form, we arrive at 

p (1 - k p 2 ) f" + [2(L + 1) - k(2L + 3)p 2 ] f + [1 - kL(L + 2)} p f = . (15) 

4 Spherical k > case 

Let us now consider the spherical case k > 0. Before starting, we recall that the coordinates 
(s, 6, <fi) we are using cover only the 'upper' half the the sphere (where r ranges in the interval 
[0, 2^])> so t ne range of s is [0, and the range of the new variable p is [0, ^]. A quite 
similar coordinate chart (with the same relation with r) covers the other half, so at the end our 
results will cover the whole sphere. 

We will prove that the equation (15) admits two different types of solutions. 

4.1 Solutions of type I 

Assuming a K-dependent power series for / 

oo 

f = J2 fnP n = f0 + flP + f2p 2 + f3 P 3 + ... 
n=0 
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then the K-dependent recursion relation leads to the vanishing of all the odd coefficients, f\ = 
fz = h = ■ • • = 0, so that / is a series with only even powers of p and a radius of convergence 
R c given by R c = 1/ \J \ H\ (determined by the presence of the second singularity). The even 
powers dependence suggests to introduce the new variable z = p 2 so that the equation becomes 



Az (1-kz) fj z + 2[(2L + 3) - 2k(L + 2)z] f' z + [l- *L{L + 2)]/ = . 



(16) 



As we have n > 0, it is convenient to complement the previously suggested variable change z = p 2 
with a further last change t = kz, with the range [0, 1] for t. Then the equation (16) reduces to 



i(l- i )/»+[(L+^)-(L + 2)t]/ t ' + ^[l-«L(L + 2)]/ = 0, 
that is, a Gauss hypergeometric equation 

t (1 - t) fa + [c - (1 + a K + b K )t] f t - a K b K f = , 

with 



(17) 



3 1 

L + -, a K + b K = L + l, a K b K = - — [1 - kL (L + 2) 



2 4k 
and the solution regular at t = is the hypergeometric function 

oo 

f(t,K) = 2 F 1 (a K ,b K ;c;t) , 2 Fi{a K ,b K ;c;t) = 1 + ^ 

n=l 

with a K and b K given by 
1 



(a K ) n (b K )n t V 



(c)n 



n ! 



a K = 



2 L 



(L + l) ± B K 



b K = 



(L + 1) T B K 



B K = 



The equation has a singularity at i = 1 that corresponds to z = 1/k (this is, to p = 1/vk or to 
r = 7r/ (2-v/k)). If the origin r = is placed in the 'north pole' of the sphere then this singularity 
is just at the equator, which is the boundary of the domain covered by the coordinate chart s. 
The property of regularity of the solutions leads to analyze the existence of particular solutions 
well defined at this point. The polynomial solutions appear when one of the two K-dependent 
coefficients, a K or b K , coincide with zero or with a negative integer number 

a K = -n r , or b K = -n r , n r = 0,1,2,... 



Then, in this case, we have 



that can be writen as 



+ 1) = -k (2n r + L + 1) , 



k/S 2 = l/((2n r + L)(2n r + L + 2)) . 



Therefore the coefficient £ that represents the sphere analogue of the modulus of the wave vector 
of the spherical wave is restricted to one of the values £ nr i given in the discrete set 



S 2 l = k (2n r + L)(2n r + L + 2) = k [(2n r + L + l) 2 - 1] , 



(18) 
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and then the hypergeometric series 2-^1 (a K , b K , c ; kz) reduces to a polynomial of degree n r in the 
variable z. 

Coming back to the equation (12), that was written making use of the radial variable s, then 
the recurrence relation leads to the following values for the even coefficients 

= kL(L + 2) -£ 2 = k{L + 2){L + A)-£ 2 = k{L + 4)(L + 6) - £ 2 

h ~ 2{2L + 3) /0 ' H ~ 4(2L + 5) H ' h ~ ' 6(2L + 7) ~ U ' " 

The result is that /, when written as a function of s, is given by the following expression: 

, kL(L + 2)-£ 2 2 [kL{L + 2) - £ 2 ][n(L + 2){L + 4) - £ 2 ] 4 
/ ~ + 2(2L + 3) S + 8(2L + 3)(2L + 5) S 

[kL(L + 2) - £ 2 }[k(L + 2)(L + 4) - £ 2 }[k(L + 4)(L + 6) - £ 2 ] 6 
48(2L + 3)(2L + 5)(2L + 7) 

which turns out to be a polynomial of degree 2n r , which will be denoted as V^ rL - 

1. Suppose that n r = 0; then £ 2 = <Sq L = kL(L + 2) and fa = 0. In this particular case the 
function / reduces to 

-P {(s) = l. 

2. Suppose that n r = 1; then £ 2 = £f L = k(L + 2)(L + 4) and fa = 0. In this particular case, 
the function / reduces to the quadratic polynomial 

f / \ 2(L + 2) 9 

3. Suppose that n r = 2; then £ 2 = £\ L = k(L + 4)(L + 6) and fa = 0. In this particular case, 
the function / reduces to the polynomial 

Pf(8)-1- 4(L + 3) K S 2 + 4 ( L + 3 ^ L + 4 ) K 2 S 4 

^2lW-1 (2L + 3) KS + (2L + 3)(2L + 5) K S ' 

4. Suppose that n r = 3; then £ 2 = £\ L = k(L + 6)(L + 8) and /s = 0. In this particular case, 
the function / reduces to the polynomial 



-/ (^ - 1 _ 6 ( L + 4 ) _2 , 12(L + 4)(L + 5) 2 4 _ 8(L + 4)(L + 5)(L + 6) 3 
3ilSj " (2L + 3) + (2L + 3)(2L + 5) (2L + 3)(2L + 5)(2L + 7) 



5. In the general case for any n r then £ 2 = £ 2 rL = k (2n r + L){2n r + L + 2), n r = 0, 1, 2, 

and then / 2 „ r / 0, but fan r +2 = 0. The funct 
2n r in the variable s, with only even powers. 



and then fan r / 0, but fan r +2 = 0. The function / reduces to a polynomial L of degree 
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Of course, one should recall that the radial solution R, is the product of s times the polynomial 
V/ t L (s). We note that in every polynomial, the coefficient c 2 i of s 2t is proportional to k 1 so that 
the 'direct' Euclidean limit k — > of all these polynomials is just the unity, that is, V nr L —> 1 
when k — y 0. However, if we come back to the equation (12) which holds for any value of k, then 
we find that for k = the recurrence relations are precisely 

-£ 2 _ -£ 2 _ -£ 2 

f2 ~ 2(2L + 3) f ° ' /4 ~ 4(2L + 5) h ' ^ ~ 6(2L + 7) /4 ' 

showing that in this case the differential equation does not admit polynomial solutions, and that 
the relevant solutions should be given by a series involving arbitrarily high powers of s (starting 
at s L ). Of course, this is to be expected as in this Euclidean case we already knew that the 
pertinent solutions of the full radial equations are the spherical Bessel functions, jiikr), which 
after extracting a factor r L are not a polynomial in the variable r. In the next subsection we 
come back on this question. 



Proposition 2 The eigenf unctions corresponding to distinct eigenvalues of the Sturm- Liouville 
problem given by the equation 

oof" + ai/' + a 2 / = 0, 

with 

a = s(l-ks 2 ), oi = 2(L + 1) - k(2L + 3)s 2 , a 2 = [£ 2 - kL(L + 2)] s , 

together with the appropriate boundary conditions in the points si = and s 2 = are 
orthogonal in the interval [0, with respect to the weight function q = S 2 ( L + 1 ) j\J\ — ns 1 . 



Proof: This K-dependent differential equation is not self-adjoint since a' / ai, but it can 
be reduced to self-adjoint form by making use of an appropriate integrating factor so that the 
equation becomes 



p(s,K)f s fAf/UM/H). (]{)) 



S 2(L+1) 



d_ 
ds 

with 

p(s,n) = s 2(L+1 Vl -ks 2 , q(s,K) = _^Z ==j \ = £ 2 - kL(L + 2). 

V 1 — K S 2 

This is a Sturm-Liouville problem that is singular since the function p(s, k) vanish in the two 
end points Si and S2] therefore, the appropriate boundary conditions are the boundedness of the 
solutions (and their derivatives) at the singular end points. The properties of the Sturm-Liouville 
problems state that even in this singular case the eigenfunctions of the problem are orthogonal 
with respect to the function q(s, k). 

The eigenfunctions are just the polynomial solutions V^ rL , n r = 0, 1, 2, . . . previously obtained 
that satisfy the orthogonality relations 

JO 
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that can be rewritten as follows 
Jo 



n r2 , k > . 



It is important to note that this last orthogonality relation is with respect a K-dependent 
weight function that coincides with measure dpk obtained in Proposition 1 for carrying on the 
quantization. 



4.2 Solutions of type II 

Let us suppose for the function f(p) the following factorization: 

f(p) = ^Jl-np 2 g(p) . 

Then Eq. (15) becomes 



p (1 -Kp ) g + [2{L + 1) - k(2L + 5)p 2 ] g' + [1 - 3£ - kL(L + 4)] pg = . 
Assuming a power series for g, 



(20) 



9 = ^2 9npr ' 

n=0 



90 + 51 P + 92P 2 + 53P 3 + ■ • ■ 



then the /-c-dependent recursion relation leads to the vanishing of all the odd coefficients, g\ = 
9s = 55 = • • • = 0, so that g is a series with only even powers of p and a radius of convergence 



R c given by R c = 1/y | n | (determined by the presence of the second singularity). The even 
powers dependence suggests to introduce the new variable z = p 2 so that the equation becomes 



4z (1 - k z) g" zz + 2[(2L + 3) - 2 k(L + 3) z] g' z + [1 - k (3 + L(L + 4))] g = . 
The change t = kz leads to 

i (1 - t) s& + [(L + ^) - (L + 3) t] g' t + -L [1 - « (3 + L(L + 4))] 5 = , 

that is, a new Gauss hypergeometric equation 

t (1 - t) g'l t + [c - (1 + a' K + b' K )t] g' t - a' K b' K 9 = , 

so that the solution regular at t = is the hypergeometric function 

00 (V 1 (h' \ t n 
g(t, K ) = 2 F 1 (a' K X;c;t), 2 F 1 (a' K X;c;t) = l + Y / ™ n ™ n 

n=l lCjn 

with a' K and 6' K given by 

1 r,, ~, „ 1 - v / /c(/e + 1) 



(21) 



(22) 



! ' 



(L + 2) ± B K 



(L + 2) T 5 K 
14 



B K 



The polynomial solutions appear when one of the two K-dependent coefficients, a' K or b' K , coincides 
with zero or with a negative integer number 

a' K = — n r , or b' K = — n r , n r = 0, 1, 2, . . . 

In this case, we arrive to the following expression for k: 

k = K /£ 2 = l/((2n r + L + l)(2n r + L + 3)) , 

that leads to the the following expression for discrete values of the energy: 

£l rl = k [(2n r + L + l)(2n r + L + 3)] = k [(2n r + L + 2) 2 - 1] . 

The recurrence relation leads to the following recursions for the even coefficients 

_ k(L + l)(L + 3) -£ 2 _ k(L + 3)(L + 5) -£ 2 k{L + 5)(L + 7) - £ 2 

92 ~ 2{2LT3) 9 °> 9A ~ 4(2L + 5) 92 ' 56 " 6(2LT^ 54 ' 

so that the function g is given by 

k(L + + 3) - £ 2 2 [k(L + 1)(L + 3) - g 2 ][^(L + 3)(L + 5) - g 2 ] 4 
ff + 2(2L + 3) S + 8(2L + 3)(2L + 5) S 

[k(L + 1)(L + 3) - £ 2 }[k(L + 3)(L + 5) - £ 2 }[k(L + 5)(L + 7) - g 2 ] 6 
+ 48(2L + 3)(2L + 5)(2L + 7) s +... 

1. Suppose that n r = 0; then £ 2 = £q L = k(L + 1)(L + 3) and 52 = 0. In this particular case, 
the function g reduces to 

V$ = l. 

2. Suppose that n r = 1; then £ 2 = £ 2 L = + 3)(L + 5) and 54 = 0. In this particular case, 
the function g reduces to the polynomial 

V? L = l- 2 ^±^KS 2 . 



(2L + 3) 



3. Suppose that n r = 2; then £ 2 = £\ L = k(L + 5)(L + 7) and g$ = 0. In this particular case, 
the function g reduces to the polynomial 



yg _ 4(L + 4) 2 4(L + 5)(L + 4) 2 4 
2i " (2L + 3) + (2L + 3)(2L + 5) ' 



4. In the general case for any n r , then £ 2 = £ 2 L = k (2n r +L+l)(2n r +L+3), n r = 0, 1, 2, . . ., 
and then g2 Ur / 0, but fl , 2n r +2 = 0. The function g reduces to a polynomial V^ rL of degree 
2n r in the variable s, with only even powers. 

The following proposition is similar to Proposition 2. 
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Proposition 3 The eigenf unctions corresponding to distinct eigenvalues of the Sturm- Liouville 
problem given by the equation 

a-og" + a 1 g' + a 2 g = , 

with 

a = s(l-ks 2 ), ai = 2(L + 1)-k(2L + 5)s 2 , a 2 = [£ 2 - 3k - kL{L + 4)] s , 

together with the appropriate boundary conditions in the points si = and s 2 = 
orthogonal in the interval [0, respect to the weight function q g = s 2 ( L + 1 )\/i — ks 2 . 

Notice that this problem is also singular (as in Propsition 2) and that the weight function 
q g is different from the one in the type I solutions. The eigenfunctions are just the polynomial 
solutions V^ rL (s), n r = 0, 1, 2, . . . previously obtained that satisfy the orthogonality relations 

Jo 

that can be rewritten as follows 

where i?f T and i?f T are the radial functions 

R n ri L = ^ Vl-KS 2 Vl iL (s, K) , Rl 2L = S L Vl-KS 2 V^S, k) . 

4.3 Final solution for the sphere and its Euclidean limit 

By working in the coordinate chart (s, 6, (ft), which covers the upper hemisphere, we have obtained 
the following solutions for the radial function R(s) in (11) (in each case, the quantum numbers n r 
and L are independent and their ranges are n r = 0, 1, 2, . . . and L = 0, 1, 2 . . .; remark also that 
had we worked in the alternative (s, 9, 4>) chart covering the lower hemisphere, the expressions 
for the solutions would have been the same) : 

(I) Solutions of type I : The radial functions are of the form 

R = s L Pi L (s,K), 

and the value of the associated energy, with quantum numbers n r and L, is given by 
£ 2 rL = « (2n r + L)(2n r + L + 2) = k N (N + 2) , N = 2n r + L. 
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(II) Solutions of type II : The radial functions are of the form 



R = S L y/l- KS 2 V l L (s,K), 



and the value of the associated energy, with quantum numbers n r and L, is given by 



Sl rL = k (2n r + L + l)(2n r + L + 3) 



kN(N + 2), N = 2n r + 1 + L. 



Two important properties are (i) in both cases £ 2 depends only on a single total quantum 
number N, so the energy levels are degenerate with respect n r and L (each level comprises 
solutions of type I and type II), and (ii) in both cases, £ 2 is proportional to the curvature. 

We have found that in the k > case, the free particle (which can be considered as a very 
special case of a central potential) is also endowed with extra accidental degeneracy, further to 
the rotational degeneracy. We also note that the expression of N as a function of n r and L is 
the same that appears in the analogous problem for the harmonic oscillator. 

Now, we can put together all the previous results and rewrite them in terms of the usual 
geodesic coordinates (r, 6, <fr). Apart of the known singularities in the angular part, these geodesic 
coordinates cover the whole sphere, with radial singularities only at r = (the origin, or North 
pole) and at r = it/(2^/k) (the antipodal point, or South pole). Now, as in both charts s = S K (r), 
it turns out that s 2 = (1 — C k {t))/k, so that in the upper hemisphere C K {r) = y/l — ks 2 and on 
the lower one C K ( r ) = — Vl — ks 2 . By reexpressing either the polynomials L (s, k) (I below) 
or the function s/l — k s 2 l (s, k) (II below) in terms of the variable £ = C K (0), we S e ^ a 
family of polynomials denoted by Q h l(0 = QnL(C K (r)) whose degrees are 

I n = 2n r , so n is even for the polynomials Q h l{0 coming from type I solutions. 
II n = 2n r + 1, so n is odd for the polynomials Qul(0 coming from type II solutions. 

and with coefficients wich are K-independent. It is clear from the previous analysis that the 
polynomials Q n L of degree n in the variable £ = C K ( r ')) provide a solution of our problem for 
all values of the radial coordinate r, this is, in both hemispheres. The symmetry of the sphere 
around its equator, conveyed by the transformation r — > tt / 1 yfH — r, reflects itself in the even 
character of the polynomials Q n i with n even and the antisymmetry of those with n odd under 
the replacement £ — > — £; notice that both type I and type II solutions are encompassed under a 
single family of polynomials Q n L (£ ) • 

The polynomials QnZ,(£)> f° r the first low values for n = 0,1,2,3,4,5,6 and arbitrary L are 
given by 



QoUO 

Q2U0 

Qsl(0 



k 2 (-l + 2(L + 2)f) 
k 3 (-3e + 2(L + 3)£ 3 ) 



1 



17 



Q 4L (0 = U (3 - 12(L + 3)£ 2 + 4(L + 3)(L + 4)£ 4 ) 
Q 5 l(6 = (l5£-20(L + 4)£ 3 + 4(L + 4)(L + 5)£ 5 ) 

Q 6 l(6 = ^6 (-15 + 90(L + 4)£ 2 - 60(L + 4)(L + 5)£ 4 + 8(L + 4)(L + 5)(L + 6)£ 6 ) 
with the following values for the global coefficients 

h = h = 2LT3 ' k4 = h= (2L + 3)(2L + 5) ' = (2L + 3)(2L + 5)(2L + 7) ' 

and in terms of these polynomials, the wavefunction corresponding to a spherical wave on the 
sphere, depending of the quantum numbers (n, L, m) is given in geodesic spherical coordinates 
(r,M)by 

*n,L, m = (S K (r)) L Q„, L (C K (r)) Y Lm (9, 0) , n = 0, 1, 2, . . . 

Notice the analogy for the radial dependence of these solutions with the angular dependence 
in the angular coordinate 9 involving the Legendre polynomials with the variable cos# and the 
associated Legendre functions involving extra powers of sin 9. This can be expected as, after all, 
the problem we are discussing is finding eigenfunctions of the Laplacian in the sphere S* 3 , which 
leads to the hyperspherical harmonics; in fact, the polynomials Q«,l(£) are proportional to the 
Gegenbauer polynomial C^ +1 (£). 

We state together several properties of the polynomials <2 n ,L(0i some of which have been 
already mentioned, while the rest can be easily derived: 



1- Qu,l(0 is a polynomial of degree n (which is even/odd for even/odd values of n). All its 
roots are real and are contained in the interval [—1,1] and as they come in pairs A, — A, 
there are precisely [n/2] roots in the interval [0, 1] and an extra root at £ = when n is 
odd (Here, [n/2] denotes the integer part of n/2). 

2. At the North pole, where £ = 1 and for all n, L, we have <2 n ,L(l) = 1; as the polynomial 
is even/odd, this implies that at the South pole, where £ = — 1 we have Q nr ,L(— 1) = ±1 
according as n is even / odd, and this for all L. 

3. At the equator, £ = 0, for n even, the value of the polynomial Q n ,z,(0) is always different 
from zero (this value is positive for [n/2] even and negative for [n/2] odd, and the absolute 
value of Q h ,l(0) decreases when n or L grows). For n odd the value of the polynomial 
Qu,l{0) is always equal to zero 



The orthogonality relations previously found, when written in terms of the radial variable r 
and of the polynomials Q n ,L(C K (r)), (with n either even or odd) become 



■K/yTR, 



S K (r) L Q ni>L (C K (r)) S K (r) L Q n2 , L (C K (r)) S K (r) 2 dr = , m ^ n 2 , k>0. 



which is a bit more general than the previous expressions, as it includes solutions of type I and 
II altogether. 
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A few polynomials with n = 4,5, 12, 13 and several values of L are displayed in the figures 1 
to 4. 

To sum up, the wave functions of the free particle in the three-dimensional sphere 5 3 , which 
have a separated expression in the geodesic spherical coordinates (r, 9,4>), depend on the quantum 
numbers (n, L, m) and are given (up to a multiplicative constant) by 

(S K (r)) L Q n , L (C K (r)) Y Lm (8, <p) , n = 0, 1, 2, . . . , 

with energies given by 

E n , L = (J^) K{n + L)(n + L + 2), n = 0,1,2,..., L = 0,1,2,... 

that are proportional to the curvature k and depend only of a total quantum number iV given by 

N = n + L = 0, 1, 2, Further to the rotational degeneracy, contained in a degeneracy factor 

2L + 1, the levels in this problem have an extra degeneracy wich leads to a total degeneracy 
(N + 1) 2 for the level with total quantum number N (e.g., for TV = 4 we have 4 = 4 + = 3 + 1 = 
2 + 2 = 1 + 3 = + 4, with total degeneracy 1 + 3 + 5 + 7 + 9 = 25, etc.) 

Another question to be discussed here refers to the k — > limit of the sphere results, where 
the known Euclidean results should be recovered. We recall briefly these Euclidean results: 
The wave functions of the free particle in the three-dimensional Euclidean E 3 , which have a 
separated expression in the geodesic spherical coordinates (r, 6, <j)) are the spherical waves that 
depend on a continuous positive label k and two quantum numbers (L, m), and are given (up to 
a multiplicative constant) together with their energies by 



^k;L,m OC j^kr) Y Lm (6, 0), 




where jL,{kr) are the spherical Bessel functions. 

The limit from the sphere to the Euclidean space E 3 , requires some care. In the polynomials 
Qn,L(,C K (,^))j whose coefficients are independent of k, all the dependence on the curvature is 
contained in the argument C K {f) which in the limit n — > goes to 1 for any value of r, so in the 
simple k — y limit the polynomial reduces to a constant. In the same 'naive' limit, with k — > 
and fixed quantum numbers, the energy would go to cero. Physically, what one should obtain 
is a continuous spectrum for the energy, bounded from below by the value E = 0; this cannot 
be done simply by means of k — > 0, but would require to maintain the value of £ in (18), what 
implies that n — > oo should be made at the same time as k — > (the quantum number L is kept 
fixed in this proccess). In other words, the simple limit k — > would make the quantum label 
n to disappear, without leaving any trace, but in the true limit, there is another quantum label 
k, which should appear. This mean that the limit should be made enforcing the constancy of 
k [(n + L + l) 2 - 1] throughout. 

When this is duly taken into account, everything fits. Of course, this is a consequence of the 
relation between the radial equation for the sphere (13) which is a kind of deformed spherical 
Bessel equation and its Euclidean limit (14) which is precisely the spherical Bessel equation. It 
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is interesting to look at this limit, where the complete Euclidean space appears as the limit of 
the upper half-sphere, with the sphere equator going to the Euclidean infinity (this fits with the 
use of the chart (s, 9, <p) which in the limit covers the complete Euclidean space). Consequently, 
the complete radial behaviour of the Euclidean radial function comes from the 'upper half, from 
£ = 1 to £ = of the full interval [—1, 1] for (C K (j)), with the equator (with r = tt/{2^/k) going 
to the Euclidean infinity as k — > 0. This corresponds to the limit 

lim (S K {r)) L Q n , L (C K (r)) oc jL(kr) (23) 

k^O; n— >oo 

whenever the control variables k (continuous and positive) and n (discrete) are linked through- 
out the limit by the condition k (n + L)(n + L + 2) = k 2 . The precise proportionality factor 
required to turn (23) into an equality comes from matching the lowest coefficients: for small r, 
(S K (r)) L Q nriL (C K (r)) « r L , while for small x, j L (x) « x L /{l ■ 3 ■ 5 ■ -(2L + 1)). 

In Figures 5 and 6 we display this fact, approaching the spherical Bessel functions quite closely 
with only a few members of the limit sequence indicated above. The concordance is good enough 
with n moderately low, which allows us to avoid the computational instabilities — catastrophic 
cancellation due to the large absolute values and alternating signs of the coefficient for the 
polynomials — which arises for higher values of n. 



5 Hyperbolic n < case 



The search of the analogues of the spherical waves in the hyperbolic space can be carried out by 
using the same approach. However, there are some relevant differences. Thus we restrict here to 
a sketch of the main traits, which will be discussed in full detail elsewhere. 

In this case, we have n = — \k\ < and by using the adimensional variables defined similarly 
to the sphere case, the equations (13) and (15) become 

p\l + |k| p 2 ) R" + P (2 + 3|k| p 2 ) R' + (p 2 -L(L + l))R = 0, (24) 

and 

p (1 + \k\ p 2 ) f" + [2(L + 1) + \k\(2L + 3)p 2 ] /' + [1 + \k\L(L + 2)] pf = , (25) 

that also represent a deformation of the Bessel equation and that can also be solved by using a 
power series procedure. The main difference is that now the singularity of the equation is placed 
on the p imaginary axis and the equation is well defined for all the real values p > 0. This fits 
of course to the fact that in the hyperbolic space, the range of the geodesic radial distance r is 
[0, oo) and this implies also that both s and p have the range [0, oo). The substitution z = p 2 
and then the change t = —\k\ z leads to the following hypergeometric equation: 
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*(!-*)/«+ (L + -)-(L + 2)t n-—[l + \K\L(L + 2)]f = 



2 



1 



4\k\ 



which nevertheless has a singularity at t = 1, which in this case comes from the singularity at 
the imaginary axis through the changes from p to z and to t. 



20 



This equation is a Gauss hypergeometric equation with c = L + 3/2 and a K and b K given by 



(L + 1) ± B K 



(L + 1) T 5 K 



5 K 



So the solution of the quantum free particle in a hyperbolic space can also be expressed in 
terms of a hypergeometric series that can be considered as representing a /-v-deformation of the 
Euclidean spherical waves. 

We close this section with two points that must be remarked. First, in the case of the hyperbolic 
space this hypergeometric equation does not admit polynomial solutions since the condition for 
a K or b K to coincide with zero or a negative integer number leads to < (a condition impossible 
to be satisfied). This fits with the physical expectations that the spectrum for the energy E = 
(h 2 /2m)£ 2 would be continuous. Second, for certain values of the curvature the coefficients a K 
and b K become complex; nevertheless, the point is that they are complex conjugated (that is, 
b* K = a K ) and because of this the solution (that is a hypergeometric function) continues to be a 
real function. In fact, we have obtained two properties (no polynomial solutions and existence of 
complex coefficients) that coincide with similar results obtained (but with a different formalism) 
in the previous study [1] of the two-dimensional system. 



6 Final comments and outlook 

We have studied the analogous to the Euclidean spherical waves for a quantum free particle 
on the three-dimensional spherical and hyperbolic spaces using a curvature dependent approach. 
We start from the K-dependent Killing vectors of the metric and we obtain the quantum (kinetic) 
Hamiltonian H(k) as a function of the operators Pi, and Jj. Then the ^-dependent Schrodinger 
equation (that when k — > reduces to a spherical Bessel differential equation) can be carried in 
the case through a number of changes of variable, to a hypergeometric equation. We study 

with full detail the case for positive curvature on a sphere, where — as the space is closed — 
quantum mechanics leads naturally to a discrete spectrum for the energy of the analogous of 
spherical waves. In the standard Euclidean quantum mechanics, this discretization is usually 
a consequence of the existence of a potential. In the spherical k > case, even without the 
presence of a potential, the geometry of the space produces a discrete spectrum. We also discuss 
in detail the limit from the sphere to the Euclidean space. 

We call the attention to some particular points. First, the differences between the spherical 
k > and the hyperbolic n < cases have been clearly stated; nevertheless, this is a point 
deserving to be studied with more detail. Second, in the spherical k > case, one can expect 
some specific relations to exists between K-dependent plane- waves (obtained in [1]) and the 
spherical waves obtained here (in the Euclidean case this a well known relation), and to find this 
explicitly can be considered as an open question. Third, in the k < hyperbolic case (that, in 
certain aspects, seems more similar to the Euclidean case that the spherical one) the equations 
can lead to complex values for the parameters a K and b K ; this also deserves be studied. 

We finalize with two interesting open questions: 
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(i) Finally, in the k > case we have obtained a family of orthogonal polynomials (that is 
different to the family obtained in [1]). We have seen in some previous papers that this 
K-formalism provides, as a byproduct, new families of polynomials; as an example, in [3] 
were obtained the so-called curved Hermite polynomials (CHP). We think that the family 
here obtained deserves also to be studied (that is, existence of the corresponding 'deformed' 
Rodrigues formula, generating function, or recursion relations) probably in a rather similar 
way as the one presented in [40]. 

(ii) The quantum free motion on the 3-dim spherical and hyperbolic spaces should correspond 
to the limiting case of the Hydrogem atom (see e.g. [6]) or the harmonic oscillator (see e.g. 
[41]) in curved spaces when the coupling constant in the potential vanishes. We note that 
in these two cases the problem involves two parameters (coupling constant and curvature) 
and there are therefore two different limiting procedures. A general analysis of these limits 
is a delicate problem that deserves to be studied. 
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Figure 1: Several polynomials Q n ,L(C K (f)) with n = 4 and values of L = 0, 1, 2, 3, 4, 5, 6, 8. 
The polynomial with L = is in dark grey, thick, and increasing L corresponds to lighter 
grey. The variable r in abscissas is the geodesic lenght on a sphere with curvature k — 1, 
so r = corresponds to the North Pole and r = n corresponds to the South pole. Notice 
all these polynomials have precisely 4 zeros and all are even around the equator, r = ir/2. 




Figure 2: Several polynomials Q n ,L{C K ( r )) with n = 5 and values of L = 0, 1, 2, 3, 4, 5, 6, 8. 
The polynomial with L = is in dark grey, thick. The variable r in abscissas is the 
geodesic lenght on a sphere with curvature k — 1, so r — corresponds to the North Pole 
and r = 7r corresponds to the South pole. Notice all these polynomials have precisely 5 
zeros and all are odd around the equator, r = ti/2 
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Figure 3: Several polynomials Q n ,L(C K ( r )) with n — 12 and values of L = 0, 1, 2, 3, 4, 5, 6, 8. 
Same conventions as in preceeding figures. Notice all these polynomials have precisely 12 
zeros and all are even around the equator, r = n/2; the amplitude of oscillations at the 
equatorial band decreases when L grows, and to allow this to be seen clearly the range 
displayed has been reduced. 




Figure 4: Several polynomials Q n ,L(C K ( r )) with n — 13 and values of L = 0, 1, 2, 3, 4, 5, 6, 8. 
Notice all these polynomials have precisely 13 zeros and all are odd around the equator, 
r = tt/2. 
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Figure 5: Aproaching the Euclidean radial function for L = (the spherical Bessel function 
jo{kr) for k — 10 with a sequence of the solutions for the sphere, corresponding to the 
values n = 20, 24, 32, 40 and L = 0. In the limit k — > the upper hemisphere goes to the 
full Euclidean plane and the sphere Equator goes to Euclidean infinity, so here only the 
left hand of the graphics in figures 1 to 4 (which corresponds to the upper hemisphere) is 
pertinent. 
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Figure 6: Same as Figure 5 but for L = 3, approaching the spherical Bessel function 
hikr). 
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